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This is a short guide to some uses of the zeta-function regularization procedure as a a basic mathematical tool for 
quantum field theory in curved space-time (as is the case of Nambu-Jona-Lasinio models), in quantum gravity models 
(in different dimensions), and also in cosmology, where it appears e.g. in the calculation of possible 'contributions' 
to the cosmological constant coming through manifestations of the vacuum energy density. Part of this research was 
carried out in fruitful and enjoyable collaboration with people from Tomsk State Pedagogical University. 



1. The method of zeta-function 
regularization 

Hawking introduced this method pj as a basic tool for 
the regularization of infinities in QFT in a curved space- 
time || H, [| . The idea is the following || . One could try 
to tame Quantum Gravity using the canonical approach, 
by defining an arrow of time and working on the space- 
like hypersurfaces perpendicular to it, with equal time 
commutation relations. Reasons against this: (i) there 
are many topologies of the space-time manifold that are 
not a product R x M 3 ; (ii) such non-product topologies 
are sometimes very interesting; (iii) what does it mean 
'equal time' in the presence of Heisenberg's uncertainty 
principle? 

One thus turns naturally towards the path-integral 
approach: 



< 92,fa,S 2 \gi,(f>i,S 1 >= J V[g, 



(1) 



where gj denotes the spacetime metric, <pj are matter 
fields, Sj general spacetime surfaces (Sj = Mj U dMj), 
T> a measure over all possible 'paths' leading from the 
j = 1 to the j = 2 values of the intervening magnitudes, 
and I is the action: 

I= \tc J( R - 2A )V=gd 4 x + J L m V=gd 4 x, (2) 

R being the curvature, A the cosmological constant, g 
the determinant of the metric, and L m the Lagrangian 
of the matter fields. Stationarity of / under the bound- 
ary conditions 



s 9\dM = °> n-dSg 



= 0, 



dAl 

leads to Einstein's equations: 
Rab - \gabR + A.g ah = 8irGT ab 



(3) 



(4) 
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T a b being the energy-momentum tensor of the matter 
fields, namely, 

1 SL m 



-9 Sg ab 



(•5) 



The path-integral formalism provides a way to deal 'per- 
turbatively' with QFT in curved spacetime backgrounds 
p. First, through a rotation in the complex plane one 
defines an Euclidean action: 



ii 



-I. 



(G) 



One can also easily introduce the finite temperature for- 
malism by the substitution ti — t± = if], which yields 
the partition function 



Z 



-pE n 



(J) 



If one now adheres to the principle that the Feynman 
propagator is obtained as the limit for f3 —> oo of the 
thermal propagator, we have shown, some time ago ||, 
that the usual principal-part prescription in the zeta- 
function regularization method (to be described below) 
is actually not needed any more, since it can in fact be 
replaced by this more general principle. 

Next comes the stationary phase approach (also 
called one-loop, or WKB), for calculating the path 
integral, which consists in expanding around a fixed 
background: 



9 = 9o + g, 



(8) 



what leads to the following expansion in the Euclidean 
metric: 



■%>0] = I[9o, M +h\g,4>] 



(9) 



This is most suitably expressed in terms of determinants 
(for bosonic, resp. fermionic fields) of the kind (here 
A, B are the relevant (pseudo-)differential operators in 
the corresponding Lagrangian): 

-l 



det 
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det 
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(10) 
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1.1. The zeta function of a fDO and its 
associated determinant 

1.1.1. A pseudodifferential operator ($DO) 

A pseudodifferential operator A of order m on a man- 
ifold M n is defined through its symbol a(x,£), which 
is a function belonging to the space S rn (R n x R n ) of 
C°° functions such that for any pair of multi-indexs a, [3 
there exists a constant C Qj/ 3 so that 

\d^a{x^)\ < C Q)/3 (1 + |£ir-H (11) 

The definition of A is given, in the distribution sense, 
by 

Af{x) = (2n)- n J j<**>a(x, 0/(0 d£, (12) 

where / is a smooth function, f € S [remember that 
S = {/ e C°°(R n );supJ:^<9 a ./(:E)| < oo, Va, S R"}], 
5' being the space of tempered distributions and / the 
Fourier transform of /. When a(x,0 is a polynomial 
in £ one gets a differential operator. In general, the 
order m can be complex. The symbol of a $ DO has 
the form 

a(x,£) = am^O+flm-i^,^ 

+o m _j(af,0 H . (13) 

being a fc (x,0 = b k (x) £ k . 

Pseudodifferential operators are useful tools j?], |8j ^| , 
both in mathematics and in physics. They were crucial 
for the proof of the uniqueness of the Cauchy problem 
and also for the proof of the Atiyah-Singer index for- 
mula. In quantum field theory they appear in any ana- 
lytical continuation process (as complex powers of differ- 
ential operators, like the Laplacian). And they consti- 
tute nowadays the basic starting point of any rigorous 
formulation of quantum field theory through microlo- 
calization, a concept that is considered to be the most 
important step towards the understanding of linear par- 
tial differential equations since the invention of distri- 
butions. 

1.1.2. The zeta function 

Let A a positive-definite elliptic ^ DO of positive order 
m € R, acting on the space of smooth sections of E, 
an n-dimensional vector bundle over M, a closed n- 
dimensional manifold. The zeta function Qa is defined 
as 

Ct(s)=tr A- s = VAT*, Res> HL = SQ . (14) 
3 

where so — dim M / ord A is called the abscissa of con- 
vergence of Ca(s) Under these conditions, it can be 
proven that Ca(s) has a meromorphic continuation to 
the whole complex plane C (regular at s = ) , provided 
that the principal symbol of A (that is a m (x, ) admits 



a spectral cut: Lg = {AG C; ArgA = 6,9% < 8 < 62}, 
Spec A<~)Lg =0 ( Agmon-Nirnberg condition) . The def- 
inition of Ca(s) depends on the position of the cut Lg . 
The only possible singularities of Ca(s) are simple poles 
at Sfc = (n — k)/m, k — 0, 1, 2, . . . ,n — l,n + 1, . . . . 
M. Kontsevich and S. Vishik have managed to extend 
this definition to the case when m € C (no spectral cut 
exists) fuj. 

1.1.3. The zeta determinant 

Let A a "J DO operator with a spectral decomposition: 
{ipi, , where I is some set of indices. The defi- 

nition of determinant starts by trying to make sense of 
the product Yii^i At , which can be easily transformed 
into a "sum": lnf^e/ Aj = Siei m Ai. From the defi- 
nition of the zeta function of A: Ca(s) = Yljcj > by 
taking the derivative at s = 0: C' A (0) = — SjgjlnAi, 
we arrive to the following definition of determinant of 
A: det f A = exp [-<a(0)] • 

A big amount of explicit formulas, useful for the cal- 
culation of zeta functions and determinants of opera- 
tors whose spectrum is known explicitly or implicitly 
(as roots of a spectral function that incorporates the 
equation and the boundary conditions, think e.g. of the 
roots of a Bessel function) are to be found, with all sort 
of explanations, in our references |ll], [l^, 0, ||, |[ pit - 
Concerning the cases of known spectrum, the most gen- 
eral situations we have been able to consider can be 
summarized as follows. They are for zeta functions of 
the form M: 

Cx(s) = £ [Q(n)+A(n)]- S , (15) 

_ rwd 

neN d 

where Q is a quadratic non-negative form and A a gen- 
eral affine form (they give rise to Epstein and Barnes 
zeta functions, respectively). I have also obtained ex- 
plicit results (given in terms of asymptotic series) for 
the much more involved cases where the summation in- 
dices are 'interchanged', namely: 

Ca(«) = E [Q(n)+A(n)]- S , (17) 

us) = Yl A ^y s - 

_ rrd 

2. Four-fermion models in curved 
spacetimes 

Four-fermion models |l4| , |l5[ — usually considered in the 
1/N expansion — are interesting due to the fact that 
they provide the opportunity to carry out an explicit, 
analytical study of composite bound states and dynam- 
ical chiral symmetry breaking. At the same time, these 
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theories — and specially their renormalizable 2d and 3d 
variants — exhibit specific properties which are similar 
to the basic behaviors of some realistic models of parti- 
cle physics. Moreover, this class of theories can be used 
for the description of the standard model (SM) itself, or 
of some particle physics phenomena in the SM. 

Having in mind the applications of four-fermion 
models to the early universe and, in particular, the 
chiral symmetry phase transitions that take place un- 
der the action of the external gravitational field, there 
has been some activity in the study of four-fermion 
models in curved spacetime. The effective potential 
of composite fermions in curved spacetime has been 
calculated in different dimensions and dynamical chi- 
ral symmetry breaking, fermionic mass generation and 
curvature-induced phase transitions have been investi- 
gated in full detail. However, in most of these cases only 
the linear curvature terms of the effective potential had 
been taken into account. But it turns out in practice 
that it is often necessary to consider precisely the strong 
curvature effects to dynamical symmetry breaking. In 
fact we could see that going beyond the linear-curvature 
approximation could lead to qualitatively different re- 
sults. 

With the help of the zeta-function regularization 
techniques, we investigated some 2d and 3d four-fermion 
models which were renormalizable — in the 1 /N expan- 
sion — in a maximally symmetric constant-curvature 
space (either of positive or of negative curvature). The 
renormalized effective potential was found for any value 
of the curvature and the possibility of dynamical sym- 
metry breaking in a curved spacetime was carefully ex- 
plored. Furthermore, the phase structure of the theory 
was also described in detail. 

One should go to the original papers for details [jl4| 
In particular, we calculated the effective potential 
of composite fermions in the Gross-Neveu model, in the 
spaces S* 2 and H 2 . The phase diagram in S 2 was con- 
structed and it was shown that for any value of the cou- 
pling constant there exists a curvature above which chi- 
ral symmetry is restored. For the case of H 2 , we showed 
that chiral symmetry is always broken. The asymptotic 
expansions of the effective potential were given explic- 
itly, both for small and for strong curvature. The three- 
dimensional case was then studied. We considered two 
different four-fermion models: one which exhibits a con- 
tinuous U(2) symmetry and another where we concen- 
trated ourselves on two discrete symmetries which hap- 
pened never to be simultaneously broken. We studied 
explicitly the dynamical P and Z2 symmetry breaking 
pattern in H 3 and S 3 . 

We started with the discussion of the Gross-Neveu 
model in the de Sitter space. This model, although 
rather simple in its conception, displays a quite rich 
structure, similar to that of realistic four-dimensional 
theories — as renormalizability, asymptotic freedom and 
dynamical chiral symmetry breaking. Some discus- 



sions of chiral symmetry restoration in the Gross-Neveu 
model for different external conditions (such as an elec- 
tromagnetic field, non-zero temperature or a change of 
the fermionic number density) had appeared in the past 
(the influence of kink-antikink configurations on the 
phase transitions was described too). 

The study of the Gross-Neveu model in an exter- 
nal gravitational field had been performed using the 
Schwinger method. Unfortunately, the generalization 
of the Schwinger procedure to curved spacetime is not 
free from ambiguities and this is why the previous re- 
sults included some mistake, which we managed to cor- 
rect, by using a rigorous mathematical treatment of the 
fermionic propagator in (constant curvature) spacetime. 
Our starting point was the action 



S 



d 2 x- 



(19) 



with N the number of fermions, A the coupling con- 
stant, 7 m (ie) = 7 a e^(a;), with 7° the ordinary Dirac ma- 
trix in flat space, and V M the covariant derivative. By 
introducing the auxiliary field a , this could be rewritten 



S 



dfx- 



i/'J7 M (x)V AI V - ^x^ 2 ~ 



, (20) 



with a = —jfijjip. Furthermore, we went to Euclidean 
space. 

There is no place for further details. Using the pow- 
erful zeta-function method, the renormalized effective 
potential was found for any value of the curvature and 
its asymptotic expansion was given explicitly, both for 
small and for strong curvature. The influence of gravity 
on the dynamical symmetry breaking pattern of some 
U(2) flavor-like and discrete symmetries was described 
in detail. The phase diagram in S 2 was constructed and 
it was shown that, for any value of the coupling constant, 
a curvature existed above which chiral symmetry was re- 
stored. For the case of H 2 , it was always broken. In 
three dimensions, in the case of positive curvature, S 3 , 
it was seen that curvature could induce a second-order 
phase transition. For H 3 the configuration given by the 
auxiliary fields equated to zero was not a solution of the 
gap equation. The physical relevance of the results was 
discussed. In particular, we could see explicitly that the 
effect of a negative curvature is similar to that of the 
presence of a magnetic field. 

For the two-dimensional Gross-Neveu model on S 2 
— where the chiral symmetry is a discrete one — we 
showed the possibility of chiral symmetry breaking and 
of fcrmion mass generation. Note that the curvature of 
a two-dimensional de Sitter space acts here as some ex- 
ternal parameter (like temperature) which induces the 
chiral symmetry phase transition. In this sense, and ow- 
ing to the fact that we treated curved spacetime exactly, 
the de Sitter space could not be considered to be some 
fluctuation over flat spacetime. 
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In the case of positive curvature, S 3 , we checked 
that the scenario given in the framework of the small 
curvature expansion (e.g., fluctuations over flat space) 
changes dramatically when gravity is treated exactly, 
what could certainly be useful in cosmological applica- 
tions. 

3. Uses in dilaton gravity 

Two-dimensional dilaton gravity interacting with a four- 
fermion model and scalars was investigated in a series 
of papers of our collaboration (see, e.g., Jl6|, 17, |l^| to 
mention only a few). The one- loop covariant effective 
action for 2D dilaton gravity with Majorana spinors 
(including the four-fermion interaction) was obtained, 
and the technical problems which appeared in any at- 
tempt at generalizing such calculations to the case of 
the most general four-fermion model described by Dirac 
fermions were properly discussed. A solution to these 
problems was found, based on its reduction to the Majo- 
rana spinor case. The general covariant effective action 
for 2D dilaton gravity with the four-fermion model de- 
scribed by Dirac spinors was given. The one-loop renor- 
malization of dilaton gravity with Majorana spinors was 
carried out and the specific conditions for multiplicative 
renormalizability were found. A comparison with the 
same theory but with a classical gravitational field was 
also done. 

Different approaches to the quantization of 2D dila- 
ton gravities (mainly, string inspired models) had been 
discussed already. We obtained the covariant effective 
action corresponding to a very general multiplicatively 
renormalizable model of 2D gravity with matter. Its 
action was of the following form 



5' 



-q{<s>)4> al x d x ^ a + &($) ($ a N ab i> h ) 



(21) 



It thus included a dilaton field, $ , n Majorana fermions 
if) a interacting quartically via a symmetric constant ma- 
trix N a i, , and m real scalars Xi ■ We also considered the 
much more difficult case in which the action contained 
2D Dirac fermions (notice that we chose the matter to 
interact with the dilaton via arbitrary functions). 

This action described and generalized many well- 
known dilaton models, e.g. the bosonic string effective 
action 



Z{§) = 8 



e 2 *, C($) = e- 2 *, = 4A 2 e- 2 *, 



<?($) = &($) = 0, /($) = 1. 
In the absence of matter our action for 
Z = 0, C($) = $, V($)=A$, 



(22) 



(23) 



coincided with the Jackiw-Teitelboim action. One could 
also add gauge fields to the matter sector. 

We constructed the covariant effective action of the 
theory (pl[), studied its one- loop renormalization and 
discussed some thereby connected issues. We described, 
in full detail, the calculation of the one-loop covariant 
effective action in 2D dilaton gravity with Majorana 
spinors, what was the first example of such a kind of 
calculation in two dimensions, as has been recognized 
latter by many authors in their mentions to our paper. 
The inclusion of scalars was also discussed there. We 
proceeded with the computation of the covariant effec- 
tive action of dilaton, scalars and Majorana spinors for 
quantum systems in classical spacetime. We discussed 
the technical problems which appeared in the deriva- 
tion of the covariant effective action in 2D dilaton grav- 
ity with Dirac fermions: the solution of these problems 
was found, via reduction of the system to the case of 
the theory of quantum dilaton gravity with Majorana 
spinors. The one-loop renormalization of quantum dila- 
ton gravity with Majorana spinors was discussed too. 
The conditions of multiplicative renormalizability were 
specified and some examples of multiplicatively renor- 
malizable dilaton potentials were explictly obtained, as 
starting point to discuss 2D quantum dilaton-fermion 
cosmology. There is no place here to describe the use- 
fulness of the zeta-function regularization procedure in 
all these developments and the reader is addressed to 
the original references, as quoted above (for a review 
of recent developments in this direction see |Ts[ |, and 
for related result with possible cosmological application 



4. Cosmological uses 
4.1. The zero point energy 

If H is now the Hamiltonian corresponding to a physi- 
cal, quantum system, the zero point energy is given by 



< 0|#|0 >, 



(24) 



where |0 > is the vacuum state. In general, after normal 
ordering we'll have: 



// = I n+^) X n a n ai, 



and this yields for the vacuum energy: 



he 



<0|#|0>=-^A 



(25) 



(26) 



(I won't normally keep track of the Ws and c's that 
will be set equal to 1.) The physical meaning of this 
energy was the object of a very long controversy, in- 
volving many first-rate physicists, until the late Hcin- 
drik Casimir gave the explanation (over fifty years ago) 
that is widely accepted nowadays, and that's the reason 
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why the zero-point energy is usually associated with his 
name. 

The expression above acquires a very important 
meaning as soon as one compares different settings, 
e.g., one where some sort of boundary conditions are 
imposed to the vacuum (e.g., a pair of parallel plates, 
infinitely conducting, in the vacuum corresponding to 
the electromagnetic field) with another situation where 
the boundary conditions (the plates) are absent (they 
have been sent to infinity) . The difference yields a phys- 
ically observable energy. 

In general the sums appearing here are all divergent. 
They give rise to the most primitive, but physically 
meaningful, examples of zeta function regularization one 
can think of. In fact, according to the definitions above: 

< 0|ff|0 >= \(,h(-1). (27) 

It is important to notice that the zero-point energy is 
something one always has to keep in mind when consid- 
ering any sort of quantum effect. Its contribution can 
be in some cases negligible, even by several orders of 
magnitude (as seems to be the case with sonoluminis- 
cence effects), but it can be of a few percent (as in sonic 
laser cavity effects), or even of some 10 — 30% as in the 
case of several wetting phenomena of alcali surfaces by 
Helium. Not to speak of the specifically devised experi- 
ments, where it may account for the full result. 

In the case of the calculation of the value of the cos- 
mological constant, it is immediate to see from the ex- 
pressions considered before that: 

<0|T - |0>= ^ + 2 J k5>' (28) 

n 

where V is the volume of the space manifold and the 
second term as a whole is the vacuum energy density 
corresponding to the quantum field (or fields) we are 
considering. Unless the first term (the cosmological con- 
stant), the vacuum energy density is not a constant (it 
goes as a -4 , a being a typical cosmological length). 
However, this does not prevent the mixing of the two 
contributions when one considers, e.g., 'the presently 
observed value of the cosmological constant'. What we 
have calculated is the second contribution for a scalar 
field of very low mass. 

4.2. A simple model 

Consider the space-time to be of one of the following 
types: RxT^xT', RxT^xS', which are actu- 
ally plausible models for the space-time topology. A 
(nowadays) free scalar field pervading the universe will 
satisfy 

(-□+£i? + M 2 )0 = O, (29) 

restricted by the appropriate boundary conditions (e.g., 
periodic, in the first case considered). We shall call p$, 



the contribution to pv from this field 

«=^=^(* 2 +^r, (30) 

where the sums and ^j^. are generalized ones (most 
common case: a multidimensional series together with 
a multidimensional integral) and p is the usual mass- 
dimensional parameter to render the eigenvalues adi- 
mensional (we take fi = c — 1 and shall insert the di- 
mensionfull units only at the end of the calculation). 
The mass M of the field will be here considered to be 
arbitrarily small and will be kept different from zero. 
This is nice, both for computational reasons as well as 
for physical ones, since a very tiny mass for the field can 
never be excluded. 

After going through some lengthy calculations that 
use the power of zeta-regularization, we reach the con- 
clusion that coincidence with the observational value 
for the cosmological constant is obtained for the con- 
tribution of a massless scalar field, p^ , for p large com- 
pactified dimensions and q = p + 1 small compacti- 
fied dimensions, p — 0, ... ,3 , and this for values of the 
small compactification length, b, of the order of 100 to 
1000 times the Planck length lp (what is actually a 
very reasonable conclusion, according also to other ap- 
proaches). To be noticed is the fact that full agreement 
is obtained only for cases where there is exactly one 
small compactified dimension in excess of the number 
of large compactified dimensions, p and q refer to the 
compactified dimensions only, but there may be other, 
non-compactifed dimensions (exactly 3— p in the case of 
the 'large' ones). In particular, the cases of pure spher- 
ical compactification and of mixed toroidal (for small 
magnitudes) and spherical (for big ones) compactifica- 
tion can be treated in this way and yield results in the 
same order of magnitude range. Both these cases cor- 
respond to (observational) isotropic spatial geometries. 
Also to be remarked again is the non-triviality of these 
calculations, when carried out exactly, what is apparent 
from the use of the generalized Chowla-Selberg formula. 
Simple power counting is absolutely unable to provide 
the correct order of magnitude of the results. 

The most precise fits with the observational value of 
the cosmological constant are obtained for b between 
b = 100 Z P and b = lOOOlp, with (1,2) and (2,3) com- 
pactified (large, small) dimensions, respectively. There 
is in fact no tuning of a 'free parameter' here. All them 
correspond to a marginally closed universe, in full agree- 
ment too with other completely independent analysis of 
the observational data |l], ||, ||, |1 • 
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